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Abstract. A solution to the a problem — that of finding 
a meehanism capable of converting Poynting energy flux 
to particle-borne energy flux in a pulsar wind — was pro- 
posed several years ago by Coroniti and Michel who con- 
sidered a particular prescription for magnetic rcconnection 
in a striped wind. This prescription was later shown to be 
ineffective. In this paper, we discuss the basic microphysics 
of the rcconnection process and conclude that a more rapid 
prescription is permissible. Assuming dissipation to set in 
at some distance outside the light-cylinder, we compute 
the resulting radiation signature and find that the syn- 
chrotron emission of heated particles appears periodic, in 
general showing both a pulse and an interpulse. The pre- 
dicted spacing of these agrees well with observation in the 
case of the Crab and Vela pulsars. Using parameters ap- 
propriate for the Crab pulsar — magnetization parameter 
at the light cylinder ctl = 6 x 10^, Lorentz factor F = 250 
— reasonable agreement is found with the observed to- 
tal pulsed luminosity. This suggest that the high-energy 
pulses from young pulsars originate not in the co-rotating 
magnetosphere within the light cylinder (as in all other 
models) but from the radially directed wind well outside 
it. 

Key words: Pulsars: general - pulsars: Crab - MHD - 
radiation mechanisms: non-thermal 



1. Introduction 

Recently, the non-axisymmetric "striped pulsar wind" in- 
vestigated b y [Coroniti (1990)| and [Michel (1994)| has been 
reexamined (Kirk & Lyubarsky 2001). A striped wind is a 
structure similar to a Parker spiral, which is produced by 
the outward radial advection of magnetic field lines whose 
foot points are anchored on the surface of the rotating 
neutron star. However, if the field lines originate from the 
poles of an obliquely rotating dipole, there exists a region 
around the equatorial plane where the polarity of the field 
at a fixed radius reverses at the rotation period of the star, 
producing stripes of alternating magnetic field direction. 
Magnetic energy is dissipated into particle energy if recon- 
nection occurs at the stripe boundaries. Coroniti (1990) 



proposed that this kind of dissipation proceeds at a rate 
sufficient to maintain the thickness of the current sheet 
approximately equal to the gyro radius of the heated par- 
ticles. Michel (1994) used essentially the same criterion, 
phrased in terms of the velocity of the current carriers. 
These early papers concluded that observations of the 
Crab Nebula, which indicate conversion of Poynting flux 
to particle-borne energy flux within lO^rL, (where tl is 
the radius of the light cylinder) could be explain ed by this 
mechanism. However, Lyubarsky fc Kirk (2001) point out 



that rcconnection is accompanied by acceleration of the 
wind. Both processes occur on the same timescale when 
the Coroniti/Michel prescription is used for rcconnection. 
This is also the timescale on which the plasma expands. 
As a result, time dilation in the accelerated flow has to 
be taken into account and the dissipation rate is signifi- 
cantly reduced — in the case of the Crab, conversion of 
Poynting flux to particle-borne energy flux is not achieved 
until roughly IO^^tl, which is too slow to account for the 



observations (the so-called cr-problem; see Melatos (1998) 
and references therein). 

In Sect. H we discuss the physics of the dissipation pro- 
cess and the prescriptions used to model it. We conclude 
that it may proceed much more rapidly than suggested 
by Coroniti and Michel. If rapid rcconnection in fact oc- 
curs, it will be accompanied by radiation losses (Melatos 
1998) that may give rise to an observable signature. In this 
paper, and in an accompanying letter (Kirk et al. 2002) 
we adopt this hypothesis. We assume that dissipation is 
triggered at a surface r = tq well outside the light cylin- 
der in the supersonic MHD wind of a pulsar and takes 
place rapidly, so that the plasma does not have time to 
expand significantly. Using two simple models of the elec- 
tron distribution we show (i) that for parameters thought 
typical of the winds from young pulsars (e.g., the Crab) 
the predicted emission is pulsed, with, in general, a pulse 
and interpulse that are not symmetrically spaced in phase 
and (ii) that the predicted synchrotron luminosity in the 
pulses agrees to order of magnitude with that observed 
from the high energy emission of the Crab pulsar. 

Our calculations suggest that the high energy pulses 
emitted by young pulsars originate as synchrotron emis- 
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sion in the pulsar wind, far outside the light cylinder. This 
is in marked contrast with conventional models of high- 
energy pulse production at an inner or outer gap — for 
a review see Harding (2001) — and also with magneto- 



spheric synchrotron based models Machabeli et al (2000) 



Crusius-Watzcl et al. (2001) The current sheet just be- 



yond the region of closed field lines in the magnetosphere 
has also been proposed as the source of high energy ra- 
diation from the Crab pulsar (Lyubarskyl996). But, in 
common with the others, this model relies for pulse pro- 
duction on the co-rotation with the star of a beam of ra- 
diation emitted from within the light cylinder. 

2. Dissipation in the current sheet 

An equilibrium solution of the Maxwell- Vlasov equations 
describing a magnetic reversal across a current sheet was 
found many years ago (Harris 1962) and was subsequently 
generalized to apply also to the case of relativistic parti- 
cle temperatures and drift speeds (Hoh 1968). If, in the 
special case in which the plasma is composed solely of 
positrons and electrons the magnetic field strength far 
from the sheet is B' , the central density of particles of 
each charge is N' and the sheet width is a, (all measured 
in the frame of reference in which the sheet is stationary 
and the electrons and positrons have equal and opposite 
drift speeds) then the temperature T, which is indepen- 
dent of position in the sheet, is given by 



T ^ 



B" 



(1) 



The Debyc radius is 



Ad = 



/ T 



V 47riV'e2 



1/2 



(2) 



and is, to within a factor of the order of unity, equal to the 
gyro radius of a particle of energy T in the magnetic field 
B' . The speed ±c/3 with which each component drifts is 
also independent of position in the sheet and is obtained 
from 



a 



(3) 



The magnetic field reverses across the sheet and the den- 
sity goes to zero far from the sheet — simple explicit ex- 
pressions for these quantities as functions of position are 
available. 

This solution is well-known to be unstable to the 
excitation of a variety of waves — see, for example. 



Daughton (1999). The standard scenario is that the cur- 



rent sheet first filaments because of the growth of the tear- 
ing mode, producing magnetic islands separated by neu- 



tral lines. These then lead to magnetic reconnection and 
its associated dissipation (Melrose 1986). This scenario 
has largely been developed in the context of the geomag- 
netic tail of the Earth. Particle-in-cell simulations appear 
to confirm its validity ([Zhu & Winglee 1996; Pritchett 
et al. 1996), although the 3-dimensional picture is more 
complicated (Biichner & Kuska 1999; Hesse et al. 2001). 

The key question in the application to pulsar winds 
is that of how quickly and, therefore, whereabouts in 
the wind the dissipation occurs. The linear growth rate 
7t of the tearing i nstability has been calculated in the 
relativistic case by Zelenyi fc Krasnosel'skikh (1979) For 
T ^ m(? they find for the most unstable mode, which has 
a wavelength approximately equal to the sheet thickness, 



(4) 



The magnetic field far from the sheet and the total num- 
ber of particles per unit area of sheet (« aN') suffice to 
determine the drift speed. A single additional parameter 
is then needed to fix the sheet properties. This can be cho- 
sen, for example, as the thickness a or the temperature T; 
the thinner (cooler) the sheet, the faster the tearing mode 
grows. 

In his model of the high-energy emission from pulsars, 



Lyubarsky (1996) argues that the rate of dissipation is 
controlled by the timescale associated with the growth of 
the tearing mode. Specifically, he assumes that magnetic 
flux is swept into the sheet and annihilated at a speed 
given by a'yt. This scenario is plausible, although it is not 
certain that the field dissipation, which presumably oc- 
curs at null points created during the field filamcntation, 
proceeds on the timescale of the tearing mode growth. If it 
operates, it results in strong radiation from the innermost 
parts of the sheet. The geometry of the flow in this region 
is, however, uncertain. 

Another, simpler, prescription for the speed 
of dissipation was (implicitly) adopted by both 
[Coroniti (1990)| and [Michel (1994)| and, later, also 
by Lyubarsky & Kirk (2001), In terms of the magnetic 
fleld and the number of particles per unit area in the 
sheet, the drift speed can be written 



/3 



B' 



SnaN'e 



(5) 



If the wind expands radially without any dissipation 
(a = constant), then its speed remains constant so that 
B' oc r^^ and N' oc r^^ . Equation (^) then implies that 
the drift speed f3 grows linearly with radius. At a radius 
where, according to this equation, /3 > 1, no equilibrium 
of the Harris type can exist; the Coroniti/Michel prescrip- 
tion consists in assuming that the sheet thickness grows at 
a rate just sufficient to avoid this happening. Dissipation 
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is thus assumed to maintain the current sheet very close to 
an equilibrium in which /? w 1. This enables reconnection 
to be incorporated straightforwardly into a set of MHD 
equations describing the pulsar wind. Provided the sheet 
does not radiate significantly, conservation of energy, mo- 
mentum and particle number, together with Eq. (||) in 
which P is set to unity, determine the dynamics of the 
wind. However, in this case dissipation proceeds on a time 
scale corresponding to the length r, which is the same 
scale on which adiabatic expansion of the sheet operates. 
A careful analysis of the dynamics shows that the wind ac- 
celerates and, as a result, the dissipation becomes rather 
ineffective (Lyubarsky & Kirk 2001). 

On the other hand, an instability in the current sheet 
can communicate itself to the rest of the pulsar wind at 
the speed of the fast magnetosonic wave, which, for a cold 
plasma is roughly equal to the Alfvcn speed. Recent in- 
vestigations of reconnection in a relativistic pair plasma 
(Zenitani & Hoshino 2001) indeed indicate that the rel- 
evant timescale is related to the crossing time of a mag- 
netosonic wave. In the case of a transsonic pulsar wind, 
a magnetosonic wavefront could appear almost stationary 
in the frame of the observer, so that dissipation could be 
completed within a few pulsar periods. This is the pre- 
scription we adopt in this paper. Energy is assumed to be 
released in the current sheet at a fixed radius, the posi- 
tion of which is taken as a free parameter. After energy 
release, wc assume the plasma cools slowly as the wind 
moves radially outwards. 

3. Geometry of the emission region 

Any radial, relativistic flow containing a periodic mod- 
ulation of the emissivity is likely to appear pulsed to an 
observer because of two effects. Firstly, the strong Doppler 
boosting of approaching parts of the flow means that only 
a small cone of the flow is visible, propagating within an 
angle of roughly l/F with respect to the line of sight, 
where F = (1 — /3^)^^/^ is the Lorentz factor and c/3 the 
radial speed. The spread in arrival times at the observer of 
photons emitted at the same radius but at different parts 
within this cone is roughly r/(cF^). Secondly, if emission 
is restricted to a range Ar in radius, the fact that the 
photons move only slightly faster than the flow when seen 
in the laboratory frame means that the spread in arrival 
times is roughly Ar/(cF^). For a wind modulated at the 
period P ~ 27rrL/c of a rotating star, the approximate 
condition for the observation of pulsed emission is there- 
fore 




Fig. 1. A meridional section of the current sheet of an 
oblique rotator with a radial wind. The thin lines limiting 
the extent of the sheet in latitude depict the plane of the 
magnetic equator when the magnetic axis lies in the plane 
of the section. 



two components (pulse and interpulse) which are symmet- 
rically placed if the pulsar is viewed from the equatorial 
plane, but are asymmetrical when viewed from higher lat- 
itudes. Denoting by a the angle between the rotation axis 
and the magnetic axis, and by C the angle between the ro- 
tation axis and the viewing direction, these pulses appear 
only for 90° - a < C < 90° a. 

To improve on the estimate (|^) and calculate the ex- 
pected pulse shapes we adopt the following model of the 
volume emissivity e of the wind: 

(i) The emission region is located on the current sheet 
in the super-magnetosonic MHD wind of a rotating 
neutron star with an oblique, split-monopole magnetic 
field. 

(ii) The dependence on frequency v of the volume emissiv- 
ity is a power-law: e oc v~'^ . 

(iii) The emission switches on abruptly as plasma in the 
sheet crosses the surface r = tq . Thereafter the volume 
emissivity decays as a power law in radius: e oc r~^~'. 

Outside the light cylinder, assumption (i) should be a 
good approximation, independent of the actual configura- 
tion of the magnetic field close to the pulsar surface. The 
assumption that emission occurs only at the exact posi- 
tion of the current sheet is an idealization — in reality, the 
entire flow may be modulated, making the observed pulse 
more complex. In this case, our computations provide the 
Greens function with which the pattern is to be folded 
to find the observed pulse shape. In the ultra-relativistic 
limit, a simple wind solution is available (Bogovalov 1999) 
in which the velocity is radial and constant and the mag- 
netic field is almost toroidal. The thin current sheet in 
such a wind is described by the equation 



A? 



< 



(6) 



= (3r\^ [± arccos (- 



cot a cot ^) 4- ct/rL — + 2n7r] (7) 



If this is fulfilled, the observed pulse shapes appear simi- 
lar to the pattern of modulation of the flow. In particular, 
for a flow containing the corrugated spiral current sheet 
of an oblique rotator (Fig. there should in general be nates, t is the observer- frame time and n is an integer. 



( Kirk fc Lyubarsky 2001 ), where (r,9,(l)) are polar coordi- 
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Assumption (ii) is an idealization that ensures the 
pulse shape is independent of frequency. In reality, the 
electron distribution may exhibit breaks and cut-offs, 
which will be translated into different pulse shapes at dif- 
ferent frequencies. Our assumption, however, ensures that 
all parts of the flow can contribute to the emission at a 
given frequency, making it harder to obtain sharply de- 
fined pulses. 

The dependence of the observed flux F ont, a and the 
direction of observation h is: 



/+00 />oc />90"+a />27r 

dt' dr d9 dcf) 

-oc Jro J90°-Q Ja 

r-iD^+"5 [r - ^{6, <j), t')] S {t' - t + h ■ r/c) 



(8) 



where D ~ 1/ (1 — /Jf • n) is the conventionally defined 
Doppler factor, apart from a (constant) factor of F. Choos- 
ing a = and q = 3, which is appropriate for the syn- 
chrotron radiation of a fixed number of particles that 
undergo predominantly adiabatic losses in the expanding 
plasma, we present sample pulses for various Lorentz fac- 
tors in Fig. H For the higher Lorentz factors this figure 
also displays the effect of a faster fall-off of the emissivity 
towards higher radii {q = 6), mimicking the effect of an 
accelerating wind. The asymmetry in time of the pulses is 
produced by our assumption of a sudden switching on of 
the emission and is reduced by a more rapid switching off. 
The parameters of this figure have been chosen to cor- 
respond to the Crab Pulsar. The radio emission of this 
object is complex, but is known to display 'core emission' 
in the radio precursor (Rankin (1990), which strongly sug- 
gests that the magnetic axis passes within a few degrees 
of the line of sight. This indicates that the inclination an- 
gle a of the magnetic axis to the rotation axis roughly 
equals the angle C between the rotation axis and the line 
of sight. The latter is determined by X-ray and optical ob- 
servations of the 'torus' (Aschenbach & Brinkmann 1975; 
Hester et al. 1995) to be 60°. In Fig. ^ the emission has 
been assumed to switch-on at ro = 30rL and pulses ap- 
pear for F > 5, in rough agreement with the estimate 
(^. At higher Lorentz factors, the pulses sharpen. From 
Eq. (^ it can be seen that the emission splits naturally 
into two contributions, corresponding to different signs 
of the term arccos(— cot a cot 6*). In the first two plots of 
Fig. ^ these contributions are plotted separately as dotted 
lines. Although the detailed pulse shapes will depend on 
the structure of the reconnecting sheets in the Crab wind, 
as well as the details of the switching on mechanism and 
the amount of acceleration of the wind, the separation of 
the two peaks is not sensitive to these unknowns, being 
determined by the location of the sheets. This separation 
in phase equals arccos(cot a cot C)/180°, corresponding to 
0.39 for the Crab Pulsar, in good agreement with obser- 
vation (e.g., Kanbach (1998) ). 




Pulse phase 



Pulse phase 



Fig. 2. Pulse shapes obtained by a numerical integration 
of Eq. (||) . Four values of the Lorentz factor (F) are shown, 
the remaining parameters are: a = 60°, C, = 60°, ro = 
30rL, a = and q = 3 (see text). The total emission is 
shown as a solid line. For F = 2 and 4, the contributions of 
the current sheets of different polarity are shown as dotted 
lines. For F = 10 and 50, the dotted line shows the effect 
of a more rapid fading of the emissivity to larger radius 
(q = 6). 

In the case of the Vela pulsar, conflicting interpreta- 
tions of the polarization sweep and the X-ray morphology 
have yielded a « 71°, C « 65° in one case (Radhakrishnan 
& Deshpande 2001) and a ^ ( ~ 53° in the other (Helfand 
et al. 2001). These imply separations of the gamma-ray 
(sub-)pulses of 0.45 and 0.31 of a period respectively, so 
that the observed spacing of 0.4 marginally favors the for- 
mer interpretation. 

4. Luminosity estimates 

Electrons heated in the current sheets will enter the sur- 
rounding magnetic field and emit synchrotron radiation. 
A rough estimate of the power radiated can be found by 
assuming a monoenergetic electron population of Lorentz 
factor 7 (3> 1) and requiring pressure equilibrium between 
the hot plasma and the surrounding magnetic field. Tak- 
ing the single particle synchrotron emissivity to be a delta 
function in frequency, the luminosity at frequency v into 
a solid angle is given by 



dL 
diydn 



/ drr^Ar^p^.p J(i. 

Jro 



1^0 ) 



(9) 



J. G. Kirk et al.: A wind model for high energy pulses 



5 



(Skjasraasen & Kirk 2001), where A (< 1) is the fraction 
of a wavelength of the wind pattern occupied by hot elec- 
trons of proper density nj^, Pg.p. = (4/3)c(TT7^(.B'^/87r) is 
the synchrotron power emitted by a single electron, and 
vq is the characteristic frequency of emission, measured 
in the observer frame (ct is the Thomson cross section). 
The condition of pressure equilibrium gives n'^jmc^ /3 = 
B'^/Stt, where B' is the magnetic field in the co- moving 
frame and m the electron mass. Well outside the light 
cylinder, but before energy release begins; i.e., at tl < r < 
To, the density and magnetic field scale as and re- 
spectively, and can be characterizes by the values tt-l and 
J5l obtained by extrapolating back to r = tl. Provided 
the flow does not accelerate significantly whilst emitting 
synchrotron radiation, these same scalings, which imply 7 
independent of r, can be used in the integration in Eq. (^). 
This is convenient for a rough estimate of the luminosity, 
but is sensitive to the assumption that the wind does not 
accelerate upon reconnection. 

If we identify the ratio of Poynting flux to kinetic 
energy (at the light cyhnder) ctl = B^/ {ATrTnLmc^), 
define an average fraction of hot particles A = 
ro drA/r^ < 1, and use the standard relation B^ = 
2.9 X io8pi/2p-5/2 ^jjg luminosity becomes 



dL 

dn = 



(10) 



where dL/df2 is the synchrotron luminosity per steradian 
divided by the total wind luminosity per steradian. Simi- 
larly, the characteristic frequency is 



/ii/o(r) = 2.6 X 10-*^ [ — j a2pl/2p-5/2 ^y^^yy (i^) 

The spectrum produced by this simple model with mo- 
noenergetic electrons extends to a maximum frequency 
given by Eq. ( pi] ) with r ~ vq and is flat: dL/diydQ oc i/". 
However, this aspect is sensitive to the way in which the 
plasma expands and accelerates after reconnection. If, in 
a more realistic scenario, the electron density is a power 
law iV(7) oc 7~^, with p > 2, then most of the energy 
resides in the particles with lowest 7 and the above es- 
timate (^l|) refers to the photons emitted by these parti- 
cles. The emission at frequencies v < voifo) is dominated 
by the same low energy particles radiating at larger ra- 
dius in a weaker magnetic field, and is determined by the 
expansion and acceleration of the hot plasma. In the ab- 
sence of acceleration, the flat spectrum derived above pre- 
vails. However, the intrinsic spectrum, which is a power 
law dL/di^dfl oc v~°', with a > 0.5, is revealed at frequen- 
cies ly > h'o^ro), where particles radiating at rg dominate. 
In the case of the Crab pulsar, the optical to X-ray emis- 
sion is consistent with a flat spectrum (Shearer & Golden 
2001), which steepens to a > 1 at a photon energy of ap- 
proximately IMeV (Kanbach 1998; Kuiper e t al. 2 001). 
Equation (|l|) therefore implies a ~ 1.1 x lO^i/ro/rL. Ob- 
servations suggest dL/dH, ~ 10^^ for this object, which, 



from Eq. (pO[), is consistent with a mildly supersonic flow, 
r « 250, that starts to radiate at rg = SOrL. This is con- 
sistent with the lower limit of 5rL found by (Bogovalov & 
Aharonian 2000) for the radius at which the conversion of 
Poynting flux to particle-born flux may occur. 

5. Discussion 

Current models of the high energy emission from rota- 
tion powered pulsars fall loosely into two groups: "po- 
lar cap" and "outer gap" models (for a recent review see 



Harding (2001) ). Each of these places the emission region 



in the corotating magnetosphere of the star. The predicted 
pulse shapes depend sensitively on the uncertain geometry 
of the magnetic field in this region, which is usually as- 
sumed dipolar. In contrast, in the model presented above, 
pulses are emitted well outside the co-rotation region. Pro- 
vided the pulsar drives a supersonic, MHD wind, the field 
geometry approaches a simple asymptotic solution (Bogo- 
valov 1999), which determines the basic properties of the 
pulses. 
In 



of 



papers. 



Pacini (1971) 



Pacini fc Salvati (1983), Pacini fc Salvati (1987) sug- 



gested that particles close to the light cylinder emit 
synchrotron radiation and are responsible for the optical 
pulses of the Crab and other pulsars — a theory which 
appears to be in reasonable agreement with observation 
(Shearer & Golden 2001). The model presented here is 
similar in that the radiation mechanism is synchrotron 
emission and the scaling of Eq. ( pO| ) is close to that origi- 
nally given by Pacini (1971), It goes further, however, by 



specifying for the emission region a precise geometrical 
structure and location. 

It is known that the dilution of the plasma in the wind 
must lead to non-ideal MHD behavior (Usov 1975, Michel 
1982) and so can be responsible for triggering the emission. 
The position at which this occurs depends on the initial 
concentration of charges in the current sheet and is taken 
as a free parameter. The major uncertainty in the model 
is the speed with which dissipation proceeds. Our compu- 
tations implicitly assume rapid dissipation (over a scale 
small compared to the radius). In this case, it is possible 
to make a rough estimate of the synchrotron luminosity 
which is in agreement with observations. A more detailed 
comparison with the wealth of relevant observations of 
optical. X-ray and gamma-ray pulsars must await further 
modeling. It is also conceivable that at least part of the 
radio emission could be produced in the wind region, as 
is suggested by the similarity of the pulse profiles at all 
frequencies seen in the Crab pulsar. However, until a can- 
didate coherent mechanism can be identified, this connec- 
tion remains very speculative. Nevertheless, the wind sce- 
nario is a viable alternative to current theories of gamma- 
ray. X-ray and optical pulses from rotation-driven pulsars 
and has the advantage of making specific predictions of 
the pulse properties. 
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